Perturbations of Keplerian Orbits in Stationary Spherically Symmetric Spacetimes 
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We study spherically symmetric perturbations determined by alternative theories of gravity to the 
gravitational field of a central mass in General Relativity. In particular, we focus on perturbations 
in the form of power laws, and calculate their effect on the secular variations of the orbital elements 
of a Keplerian orbit. We show that, to lowest approximation order, only the argument of pericentre 
and mean anomaly undergo secular variations; furthermore, we calculate the variation of the orbital 
period. We give analytic expressions for these variations which can be used to constrain the impact 
of alternative theories of gravity. 



I. INTRODUCTION 



o . 

In recent years, there has been much interest in studying theories of gravity alternative to General Relativity (GR): 
these theories have been motivated by the current observations, which seem to question the GR model of gravitational 
interactions on large scales, such as the galactic and cosmological ones. Just to mention some examples, one can 
consider f(R) theories of gravity [H-[l] or Modified Newtonian Dynamics (MOND) Q. Other theories alternative to 
GR are motivated by the attempts of giving a description of the gravitational interactions in a quantum framework: 
the recently proposed Hofava-Lifshitz gravity 0-0] is one example. 
' In order to get a deeper insight into these theories, it is important to test their predictions in a suitable weak-field 
and slow-motion limit: to this end, for instance, central mass solutions which generalize the one obtained in GR (the 
i Schwarzschild solution) are investigated and their predictions compared to the data available from astronomical and 

astrophysical observations, in the Solar System and beyond. For these purposes, a general theoretical framework 
that can deal with classes of metric theories of gravity has been developed: the parametrized post- Newtonian (PPN) 
formalism. In the latter, the weak-field and slow-motion limit of these theories is studied in terms of suitable param- 
eters, in such a way that experiments may fix their values: as a matter of fact, the current best estimates of the PPN 
parameters are in agreement with GR Q. As a consequence, it is reasonable to believe that the effects of theories 
■ alternative to GR are small, so that they can be dealt with as perturbations of the GR background. 

In this paper we want to suggest a simple approach that, without requiring the complex and comprehensive PPN 
framework, allows to evaluate the predictions of alternative theories of gravity on the Keplerian orbit of a test particle, 
which can be thought of as a simplified model of the dynamics of celestial bodies in planetary systems. In particular, 
we consider a general stationary and spherically symmetric (SSS) spacetime metric, that can be thought of a solution 
of the field equations of a generic alternative theory of gravity describing the gravitational field around a point-like 
central mass, and we work out the perturbations of the Keplerian orbital elements. We focus on perturbations in 
form of power laws, which have been recently considered in the literature [sj]: these perturbations are interesting 
because they reproduce some known solutions of the field equations of alternative theories of gravity; moreover 
5^ , arbitrary spherically symmetric perturbations can be written in terms of power series, so that our results can be used 
quite generally. In particular, by using the Gauss perturbation scheme, we obtain, to lowest approximation order, 
the expressions of the secular variations of the orbital elements and the orbital period in terms of hypergeometric 
functions. These expressions can be used to place bounds on the parameters of alternative theories of gravity. 

The paper is organized as follows: in Section [TT] we write the perturbing acceleration in a generic SSS spacetime 
and obtain its expression to lowest approximation order, while the Gauss perturbation equations are briefly reviewed 
in Section [IIII In Section HVl we focus on perturbations in the form of power laws, while conclusions are in Section fVl 
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II. THE PERTURBING ACCELERATION 



We suppose that the gravitational field of a point mass in a generic alternative theory of gravity is described by 
the SSS spacetime metric 1 

ds 2 = (1 + 4>{r)) dt 2 - (1 + iP(r)) (dr 2 + rW) (1) 

where dfl 2 = dQ 2 + sin 2 8d<p 2 , and <fi(r), ip( r ) are functions depending on the mass M of the source and, possibly, on 
other parameters of the theory. The metric ([l} is written in isotropic polar coordinates, such that the spatial part of 
the metric is proportional to the flat spacetime metric dr 2 + r 2 dfl 2 = dx 2 + dy 2 + dz 2 . Consequently, it is possible to 
write the metric ([T]) in Cartesian coordinates in the form 

ds 2 = (1 + 4>{r)) dt 2 - (1 + VO)) (dx 2 + dy 2 + dz 2 ) (2) 



where r = y 1 ' x 2 + y 2 + z 2 . 

Since the effects of the generic alternative theory of gravity are expected to be small, we suppose that they can be 
considered as perturbations of the known GR solution, i.e. the Schwarzschild spacetime. This amounts to saying that 
4>(r) and ip( r ) m some limit have to approach their GR values. In other words we suppose that in suitable limit, the 
metric takes the form 

4>{r) = <f>Gn(r) + Mr) (3) 
ip(r) = ipGR(r) +ipA(r) (4) 

where the GR values are given by 4>gr{ t ) = -2M/r + 2M 2 /r 2 , ipGR,( r ) — 2M/r (see e.g. [l(|) and the perturbations 
4>a{t), iPa(t) due to the alternative gravity model are such that (j)A(r) <C 4 } gr{i'), ipA(r) <§C V'GflX 7 ")- 

In order to calculate the perturbations of the orbital elements, we must calculate the perturbing acceleration. To 
this end, first we assume that in the given theory, the matter is minimally and universally coupled, so that test 
particles follow geodesies of the metric (|2|) (or, equivalently (Q}). Then, we consider the (post-Newtonian) equation 
of motion of a test particle (see [ll|) 

x l = ~h 00 ,i - ^h lk h 00: k + h 00tk x k x l + (h lk , m - ^fcm,i^ x k x m (5) 

where "dot" stands for derivative with respect to the coordinate time. Since in our notation it is — <j)(r) and 
hij = %jj(r)8ij, we can write the perturbing acceleration W in the form 

W = — {$(?•) [1 + ip A {r)] + *(r> 2 } x + [$(r) + *(r)] (x • v) v (6) 
where we set x = (x, y, z), v = (x, y, i), x = x/|x|, and 

dr dr 

We aim at investigating the lowest order effects on planetary motion of <pA( r ) and ijjA( r )'- to this end, it is sufficient 
to apply the Gauss perturbation scheme to a Keplerian ellipse. 

Moreover, since we are interested in the lowest order effects, we may also neglect the non linear terms (i.e. non 
linear perturbations with respect to flat spacetime). To this end, we start by noticing that to Newtonian order, it is 
v 2 ~ <t>GR( r ) — M/r. As a consequence, in © we may neglect the terms proportional to v 2 and to (x • v) v (which 
is also proportional to the orbital eccentricity) and also the term $(r)^(r). In summary, in the weak-field and 
slow-motion limit the perturbing acceleration that we are going to consider is purely radial and is given by 

W = --$(r)x = W r it (8) 

We point out that even though we started from the study of a SSS spacetime with the aim of setting bounds on 
alternative theories of gravity, what follows can be applied to generic radial perturbations (in particular, to those in 
form of a power law) of the orbital elements of a Keplerian motion. 



1 If not otherwise stated, we use units such that c = G = 1; arrowed bold face letters like x refer to spatial vectors while Latin indices 
refer to spatial components. 
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III. PERTURBATION EQUATIONS 



We start from the expression (JS) of the acceleration and apply the Gauss perturbation scheme to derive the secular 
variations. Because of spherical symmetry, the motion of test particles is confined to a plane and, in this plane, the 
unperturbed Keplerian ellipse is 



a (l 



1 + e cos / 



(9) 



where a is the semimajor axis, e the eccentricity, / the true anomaly describing the particle angular distance from the 
pericentre. For purely radial perturbations, the Gauss equations for the variations of the Keplerian orbital elements 
read [3 
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where M is mean anomaly, / is the orbital inclination, is the ascending node, uj is the argument of pericentre, 
n = U M/a 3 is the Keplerian mean motion. In a Keplerian orbit the mean anomaly is a linear function of time defined 
by M. = n(t — to), where to is the time of a passage through pericentre. The orbital period Pb is related to the mean 
motion by n — 2n/Pb. 

In order to obtain the secular effects we must evaluate the Gauss equations onto the unperturbed Keplerian ellipse 
(j9]), and then we must average them over one orbital period of the test particle. From (fl0|) and (fTTj) . and taking 
into account ([5]), we see that the secular variations of the semimajor axis and eccentricity are null, because when 
averaging them the arguments of the integrals are odd functions. Hence, we obtain that in SSS spacetimes, to 
lowest approximation order (i.e. when the perturbations are purely radial) semimajor axis, eccentricity, beside the 
inclination and the node (which are only affected by normal i.e. out-of the-plane perturbations), are not affected 
by secular variations. We point out that, once that the secular variations of the argument of the pericentre < uj > 
and mean anomaly < M. > have been determined, it is possible to obtain the corresponding variation of the mean 
longitude < A >=< uj > + < M > [Ill- 
Starting from the equation describing the variation of the mean anomaly, it is possible to evaluate the perturbation 
of the orbital period. In fact (see [13)) on using 



df = ^ 2 sf\ - e 2 dM 



from (fTS"|) it is possible to write 
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Since we suppose that the perturbations are small, we may write 

1 
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And then, by integrating over one revolution along the unperturbed orbit, we obtain 



P~P b + P A 



(19) 
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where 



is the unperturbed period, and 

Pa 



P h = 
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-W r cos/ 



(20) 



(21) 



is the variation of the orbital period due to the perturbation. 

The above equations (jT4"]) . (IT51) and (|2~T1) allow to calculate the variations of the argument of pericentre, mean 
anomaly and orbital period. This can be accomplished at least by numerical methods for arbitrary perturbations, 
however, as we are going to show in next Sextion, when the perturbation is in the form of a power law, it is possible 
to obtain analytic expressions. 



IV. POWER LAWS PERTURBATIONS 



After having discussed in the previous Section the expressions for the secular perturbations of the orbital elements 
for a generic purely radial perturbation, here we focus on the particular case of perturbations in the form of a power 
law. As we are going to show, in this case we can obtain analytic expressions for the secular variations. We point 
out that an arbitrary spherical symmetric perturbation that is expressed by an analytic function can be expanded in 
power series up to the required approximation level: hence, by knowing the contribution of each term of the series, 
it is possible to evaluate the whole effect of the perturbation, within the required accuracy. Eventually, even though 
our approach is motivated by the study of the effects of alternative theories of gravity, the results that we are going 
to obtain are quite general, and apply to radial perturbation of Keplerian orbits by means of arbitrary power laws. 

In particular, we consider two kinds of perturbation that we write as follows. Given a constant a, which is a 
parameter deriving from the gravity model alternative to GR, the perturbations that we focus on are in the form: 
(i) 4>a{t) = tw where the integer N is such that N < —1, or differently speaking 4>a(i~) — ^wji \N\ > 1; (ii) 
4>a{t) = aT " N where the integer N is such that N > 1. The perturbations of the first kind are asymptotically 
flat, while the second ones are not: in the latter case, we assume that there exists a range < r < f for which 
4>a(t) <C (/)Gij(r), and our results are valid within this range. We notice that, for a logarithmic perturbation in the 
form 4>A{r) = /3 log(r/ro), it is W r — —57, which can be dealt with as in the case (i) above. 

The expressions of the secular variations that we give below can be used to place constraints on a by means of a 
comparison with the available data. 



A. Perturbation of the argument of pericentre 

For perturbations in the form <pA (r) — -^rj, \N\ > 1, we start from eq. (|14[) . and we average it over one orbital 
period by taking into account the expression of the unperturbed orbit ((9|) and the relation [l4| 

fl - e 2 ) 3 / 2 
nil + ecosjY 

On using (IA3|) . we then obtain 

l7ra|JV|(|JV|-l)(l-e 2 ) 1HAr| / |iV| 3 |iVl 2 



< ^ >= 2 nW+m F I 1 - V' 2 - V' 2 < G ) ^ 

where F is the hypergeometric function (see Appendix . In particular, this result is in agreement with the corre- 
spondent expression found in Q , where the effects of a central force were considered. 

On the other hand, for perturbations in the form 4>A{f) — arN with N > 1, we start again from (|14j) but, in order 
to average it over one orbital, it is useful to introduce the expression of the unperturbed orbit 



r = a(l — ecosE) 



(24) 
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in terms of the eccentric anomaly E, which is is related to the true anomaly / by the following relations: 



On using the relation \U 



n 

and the integrals (IA3|) and (| A4|l . we obtain 



. cosE'-e . , \J 1 - e 2 sin E 

cos/ = - -, sm/= — — . 25) 

1 — e cos E 1 - e cos E 



dt = - (1 - ecos E)dE (26) 
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(27) 



or, equivalently 



l TraAHiV + lh^g / TV 1 TV 
< w >--2 F [ 1 --2'2--2' 2 ' e ' (28) 



In particular, eq. (|28[) is in agreement with the correspondent expression found in [9j. 

These general expressions can be used to reproduce known results. For instance, on setting TV = 1 in eqs (|28[) we 
obtain the variation of pericentre due to a constant perturbation acceleration 



71 IX I 

<uj>= —y/i-e? 29 



This result is in agreement with Sanders (also reported by @), who, working in MOND framework, considered the 
case of the constant acceleration to put constraints on the effects of anomalous acceleration by analyzing Mercury's 
advance of perihelion. 

The perturbation in the form 4>A{r) — ar 2 coincides with the effect of a cosmological constant A [16], [17] in GR, 
while in the case of f(R) gravity, it describes the vacuum solution in Palatini formalism [3, In this case, by 
setting N — 2 in eq. (f28| we get 

3na 



< u >= - — VI -e 2 (30) 
On setting a = — iA and using the relation n 2 a 3 = M which holds for the unperturbed orbit, eq. (|30j) becomes 



3 



<w>=^VW> (31) 

in agreement with [16] and [19j . 

On setting \N\ = 4 in (j2"3")l . we obtain the perturbation determined by a vacuum solution of Hofava-Lifshitz gravity 
[20L [2ll | ; by approximating the hypergeometric function, the variation of the argument of pericentre becomes 



6wa (1 + \e 2 ) 
W (1-e 2 ) 3 



< w >~ —± i-£ 32 



in agreement with [20j. 

Eventually, we consider a logarithmic perturbation, in the form ^(r) = /31og(r/ro). From eq. (|23j) . by approxi- 
mation up to e 4 we obtain 



2 n 2 a 2 V '2' ' 7 2 

in agreement with 0. 
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B. Perturbation of the mean anomaly 

We may proceed as in the previous section to calculate the secular variation of the mean anomaly. For a perturbation 
in the form <Pa(i") — 71^71 1^1 > lj we start from eq. (IT51) . and focus on the part that is proportional to W r ; after 
averaging over the unperturbed ellipse © making use of eqs. (|2"21 and (| A4|) . we get 



< M >-- 



2na\N\ (l 



> 3/2-|W| 



,a fl a+|Ar 



2 _ ~2~' 2~' ,e )-vl-e 2 <w> 



(34) 



where < v > is given by ([23]). 

Similarly, for perturbations in the form <pA{ r ) = ar N , with N > 1, averaging the part of (|I5[) which depens on 
making use of (|26|) and (|A4j) . we obtain 



< A4 >= 



2?raA 



TV 1 N 



(35) 



where < u > is given by (|28|) . 

On setting N — 2, from eq. (|35[) . by approximating the hypergeometric function, the variation of the mean anomaly 



is 
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< M >= — s-jF -l,--,l,e 2 



37ra ,„ 9 n 37ra ( 7 9 
— =- (1 -e 2 ) ~ — =- - +e 2 



(36) 



in agreement with |19l |. 



C. Perturbation of the orbital period 



Starting from eq. (1211) . we may calculate the variation of the orbital period. In particular, for a perturbation in 
the form 4>a{t) = -jw\, \N\ > 1, on evaluating eq. (|2Tj) over the unperturbed ellipse (HJ) making use of the integrals 
([A"3|l . JSU, we obtain 



Pa 
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n 3 a 2+|JV| 



(37) 



As for perturbations in the form 4>a{t) = ar N with N > 1, on evaluating eq. (|21[) over the unperturbed ellipse ([23 
making use of the relation (fT6|) and the integrals (|A3[) . (|A4|I . we obtain 
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(38) 



On setting |AT| = 4 , from eq. (|37p . by approximating the hypergeometric functions, we obtain 



Pa 



ira ( 1 — e 



2 + -e 2 



(39) 



in agreement with [2l| . 

Furthermore, on setting A" = 2 in (|38[) . by approximation of the hypergeometric functions, the perturbation of the 
orbital period turns out to be 



in agreement with (l7j . 



P^-™(7 + 3e 2 



(40) 
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V. CONCLUSIONS 



We considered a general stationary and spherically symmetric spacetime, and worked out the perturbations deter- 
mined by a generic alternative theory of gravity to the GR solution describing the gravitational field around a central 
mass. In order to evaluate the effects of such perturbations, we showed that, in the weak-field and slow-motion limit, 
to lowest approximation order, these effects can be described by a purely radial acceleration. Then, we considered the 
Keplerian orbit of a test particle, which can be thought of as a model of the dynamics of celestial bodies in planetary 
systems, and we evaluated the impact of such perturbations on the orbital elements. 

In particular, we obtained analytic expressions, in terms of hypergeometric functions, for the secular variations 
of the advance of pericentre, mean anomaly and orbital period determined by perturbations in form of power laws; 
the other orbital elements do not undergo secular variations to lowest approximation order. The expressions for the 
variation of the argument of pericentre are in agreement with the results recently obtained, pertaining to the orbital 
precession due to central force perturbations. 

Our results are quite general, since even though were motivated by the study of the effect of alternative theories 
of gravity, they can be applied to arbitrary radial power laws perturbations of the orbital elements of a Keplerian 
motion. 

Spherically symmetric perturbations in the form of a power law were obtained both in GR, for instance when the 
effect of a cosmological constant is considered, and in alternative theories such as f(R), MOND, Hofava-Lifshitz 
gravity; we have shown that our results are in agreement with those already available in the literature pertaining to 
such alternative theories. For an arbitrary perturbation it is possible to obtain the secular variations by means of 
numerical approaches, or by expanding it in power series and applying our results to each term of the series. 

The possibility of testing the predictions of alternative theories of gravity on planetary motion is important to 
verify their reliability on scales different from those typical of galactic dynamics or cosmology, where they are usually 
tested. The simple approach that we considered here allows to place bounds on the theory parameters, by evaluating 
their impact on the Keplerian dynamics and allowing a comparison with the available data. 



Appendix A: Solutions of Integrals by means of Hypergeometric Functions 

In order to evaluate the secular variations of the Keplerian elements and of the orbital period determined by the 
perturbations in the form of power laws that we have considered, it is necessary to solve the following integrals 

I N = I cos u [I + e cos u] du (Af) 
Jo 

p2-rr 

L N = / [l — ecosu] N du (A2) 
Jo 



They can be solved by subsituting cosu = z, and then using the binomial theorem for (1 + z) . The solutions are 
given in terms of hypergeometric functions: 

J w = ,JV e F^-^l-^2,e 2 ) (A3) 



where F(a, b, c, x) = 2 i*i(a 3 6, c, d) defined by (see e.g. [221 ]) 

OO / 

\- (a 

(c) n n\ 



2 F 1 (a,b,c,x) = 2_^—p, p (A5) 



n=0 

where the Pochhammer symbol (a) n is defined by 



(a + n- 1)! 
(a -1)1 

(o)o = 1 (A7) 



(a)« = (a _ 1) / (A6) 
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